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ABSTRACT 

We compute the deflection angle to order (m/ro)^ and m/rg x Kr^ for a light 
ray traveling in a flat ACDM cosmology which encounters a completely condensed 
mass region. We use a Swiss cheese model for the inhomogeneities and find that 
the most significant correction to the Einstein angle occurs not because of the 
non-linear terms but instead occurs because the condensed mass is embedded 
in a background cosmology. The Swiss cheese model predicts a decrease in the 
deflection angle of ~ 2% for weakly lensed galaxies behind the rich cluster A1689, 
and that the reduction can be as large as ~ 5% for similar rich clusters at 2; ~ 1. 
Weak lensing deflection angles caused by galaxies can likewise be reduced by as 
much as ~4%. We show that the lowest order correction in which A appears is 
proportional to m/rg x \/Atq and could cause as much as a ~0.02% increase 
in the deflection angle for light that passes through a rich cluster. The lowest 
order non- linear correction in the mass is proportional to m/ro x ^JmjrQ and can 
increase the deflection angle by ~ 0.005% for weak lensing by galaxies. 

Subject headings: gravitational lensing — cosmology: theory 



1. Introduction 



Recently iRindler fc IshakI (j2007| ) have stirred interest in the possibility of measuring the 
cosmological constant A through its effect on the deflection of light that traverses large galaxy 
clusters by asserting that A has a non-negligible effect on sm all angle bending. Several pa- 



pers have since app e ared to support the ex istence of an effect (Ilshakll2008ulshak et al.ll2008 



2OIOI : ISereno I l2008l . l2009l : ISchiickerl l2009af) although qualitativ ely disagreeing on its value 
and/or interpretation, but others ( lParkll2008l : lKhriplovichll2008l ) contest the existence of an 
effect arguing that the additional bending caused by A vanishes when measured by observers 
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moving with the Hubble flow. We purport to give the definitive answer to this question as 
well as several other related ones. When comparing observations with and without a A one 
must compare observations of two different sets of events, by two different observers, in two 
different universes. One ideally attempts to make common as many kinematic and dynamic 
properties as possible in the two gedanken experiments. To conclude whether A does or 
doesn't cause bending can easily depend on what is held in com mon and what property is 



compared in the two experiments. For example, a photon orbit in a lKottlerl (119181 ) spacetime 
(Schwarzschild with a cosmological constant) does not depend on A if static coordinates are 
used, se e Eqs. ([3|) and (flT]). On e could hence conclude that A does not effect bending. How- 



ever, as 



Rindler &: IshakI (120071 ) point out, observers, stationary relative to the Schwarzschild 
mass, will measure an angle between the photon direction and the radial direction that does 
depends on A. From this observation, one could conclude that A does effect bending. Both 
conclusions are valid, but neither answers the outstanding question, "How does A contribute 
to deflections caused by large inhomogeneities in the otherwise homogeneous background cos- 
mology?" This is just one of the questions we deflnitively answer, subject to the condition 
that the inhomogeneity is signiflcantly condensed and has no peculiar motion. 

To correctly analyze A's non-linear effect on bending we found it necessary to use exact 
solutions to Einstein's equations. These solutions reveal a somewhat surprising value for 
the lowest order non-linear correction in the deflecting mass m to the familiar Einstein 
deflection formula AGm/c^ro [see the square root term in Eq. ( !32|) ]. This correction, like 
the A correction, increases the deflection and occurs because the deflector is embedded in a 
universe which expands. By using an exact inhomogeneous cosmology the largest correction 
to AGm/c^ro is revealed not to be a non-linear term but instead is caused by the limited 
time the deflector has to influence a passing photon [see the cos^ 0i term in Eqs. (132|) and 
( 133|) ]. The limited time or equivalently the limited range of the inhomogeneity can be 
thought of as a result of shielding by the homogeneous background in which the deflecting 
mass is embedded and decreases the deflection (relative to the Einstein value). General 
relativity (GR) requires that the two gravity flelds, the homogeneous background and the 
local inhomogeneity, be appropriately matched at their bounding surface. Because the Swiss 
cheese models are the only known exact GR solutions that embed spherical inhomogeneities 
in expanding homogeneous universes, we use them. Gravitational lensing calculations in 
cosmology are usually done by superimposing a deflecting mass on top of a homogeneous 
mass density and ignoring any boundary matching. The resulting deflection angle is obtained 
by a simple summation of the Einstein expression AGm/ c^tq. The results can at best be 
accurate to flrst order in the ratio of mass m to minimum impact distance tq, and if the 
shielding predicted by boundary matching in Swiss cheese is accurate, the linear term can 
be in error by a few percent in physically realizable circumstances (see Figure 4). 
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Because our goal is to correctly present the higher order corrections and because a simple 
superposition of the masses isn't satisfactory, we resort to using exact solutions to Einstein's 
gravity (see Section 2). It is in the non-linear corrections that the cosmological constant first 
appears. What we calculate in this paper is (in a series approximation) the angle a = ^2 — 
(see Fig. 1) between two spatial directions ^2 and ^1 as seen by co-moving observers in a flat 
Robert son- Walker (RW) spacetime, where ^1 is the spatial direction of a photon before it 
encounters an inhomogeneity here described by the vacuum Kottler spacetime [see Eq. ([3])] 
and ^2 is the direction of the photon after it has emerged from the Kottler condensation. 
The dynamics of the RW metric is determined by General Relativity sourced by pressure 
free matter (often referred to as dust) and includes a cosmological constant. These are 
relatively simple Friedman-Lemaitre-Robertson- Walker (FLRW) cosmological models, see 
Eq. dl]). Because the RW cosmology used is spatially flat and non-rotating (a) the angle 
between any two co-moving spatial directions is well defined independently of when or where 
the directions are measured, and (b) the spatial direction of an undefiected photon remains 
fixed. 

In Section 2 we describe the inhomogeneous cosmology we use and in Section 3 we 
outline some details of how we compute the bending angle of a photon that encounters an 
inhomogeneity. In Section 4 we discuss limits on the usefulness of our results and compare the 
Einstein angle with our corrected results for deflections caused by inhomogeneities ranging 
from galaxies to rich clusters. 



2. Swiss Cheese: Locally Inhomogeneous Cosmologies 



We use a single condensation in a Swiss cheese cosmology to compute light deflec- 
tions caused by local inhomogeneities so there can be no doubt about errors introduced 
by gr avity approximations. Because Swiss che ese is an exact solution to Einstein's equa- 
tions (lEinstein fc Straus! Il945l : ISchiickingl |l95J) its use insures the accuracy of the super- 
posed gravity field and automatically takes into account the finite range of the mass den- 
sity pe rturbation as well as observer aberration. The model was first used by iKantowski 
(Il969al ) without the cosmological constant to settle the dispute as to whether or not inho 



mogeneities effected mean luminositi e s. At that time, the validi t y of prediction s (IZel'dovich 



1964J : iDashevskii fc Zel'dovichl Il965l : IPashevskii fc Slvshl Il966l : iBertottil Il966l : iGunnI 119671 ) 
obtained using approximate GR solutions, which indicated that mass concentrations caused 
the average distance-redshift relation to differ from the pure homogeneous value, were being 
questioned. Even though the results are on occasion dou bted by some, the Swiss cheese 



model gave the definitive answer, there is an effect, see IKantowski. Vaughan. fc Branch 
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( 119951 ) ■ This model again comes to the rescue by clearly demonstrating the extent to which 
the cosmological constant A influences the small angle bending of a photon that passes a 
single mass concentration, see Eq. fl5^ . Even though others have computed bending angles 
that depended on A, until now, questions abound as to their usefulness and/or accuracy in 
a cosmological setting. We have succeeded in giving a rigorously derived expression for this 
deflection. 

The Swiss cheese model simultaneously accounts for the finite size of the deflector's in- 
fluence, the motion of the cosmic observers, and the non-linear effects of gravity. The basic 
idea of Swiss cheese is to remove non-overlapping co-moving spheres of homogeneous dust 
from one of the homogeneous FLRW cosmologies and replace them with gravity fields repre- 
senting appropriately condensed spherical mass distributions. If the cosmology is without a 
cosmological constant the simplest replacement s are Schwarzschild me trics and if there is a A 
the simplest replacements are Kottler metrics (jPyer fc Roederlll974l ). These condensations 
are the most extreme for Swiss cheese. An infi nity of less extreme models can be construc ted 



by using the Lemaitre-Tolman-Bondi metrics (lLemaitrelll933t lTolmanlll934l : iBondil 119471 ) to 



represent spherically symmetric dust concentrations (jKantowski 1969b ). In all Swiss cheese 



models the metric that is used to fill a dust condensation must match first and second fun- 
damental forms on the boundary. In the case of Schwarzschild the metric's mass is fixed 
by the dust's density and the size of the condensed hole, and in the case of Kottler, A is 
additionally required to be the same inside as out. In this calculation we stick with the 
extreme but uniq ue condensation, the Kottler metric, to arrive at a unique deflection angle. 
Schiickerl (j2009bl ) uses this same model but because he only considers a single numerical 
example his results are difficult to compare with ours. 

The two metrics are: outside, flat (^a + Qm = 1) FLRW 

ds^ = -c^dT^ + R{Tf[dx^ + x\de^ + sin^ed(P% (1) 



with the cosmic time development given by the Friedman equation 



R 
R 




Rq 

'r 



and inside, the static Kottler metric (lKottlerlll918l ) 

ds"" = --f{r)-^c^dt^ + -iir fdr'' + r^dO^ + sin^ 6 dcj?) 
where 7(r) is defined by 

Ts Ar^ 
r 3 



7(r)^l/Wl 



(2) 



(3) 



(4) 
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Boundary matching at co-moving FLRW radius Xb constrains the Schwarzschild radius of 
the condensed mass to be ^ 

rs^Qm^iB^Xbf, (5) 



c2 



and the additional Kottler parameter A to coincide with the FLRW value, i.e., 



H, 



A^SQa^. (6) 



The Kottler and RW angular coordinates are matched at the boundary and the radius of 
the Kottler hole expands according to 

n{T) = R{T)xi, (7) 

(for some numerical examples see the Mass and columns of Table 1). As seen by a station- 
ary Kottler observer the dust boundary of the Kottler hole moves with Lorentz parameters 
/3h and given by 

76 = (8) 

-2 



= yl-Tb ^HbTb/c, 

where Hb is the time dependent Hubble parameter of the boundary. The normalized 4- 
velocity of the boundary coincides with the RW co-moving dust velocity urw at the boundary 
and is of the form 

Urw = 'lbUK + Pblb rK, (9) 
when expressed in terms of unit Kottler time and radial vectors, respectively uk and rx- 



3. The Photon's Path 

In Fig. 1 we show the spatial orbit of a shghtly deflected photon r(0) that enters and 
exits a Kottler condensation. The coordinates have been rotated to put the orbit in the 
9 = Ti/2 plane and to make it symmetric about = 7r/2 while in Kottler. The tangent to 
the photon's geodesic path is 



To 



7o 



Uk ± 




'0 ''^0 1 

- (Pk 



(10) 



where I is an angular momentum like constant; tq is the minimum r and occurs at = 7r/2; 
7o = 7(ro); and uk^ vk and 4>k are unit vectors pointing respectively in the static time. 
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radial and azimuthal Kottler directions. The actual orbit r(0) is approximated as 



r/ro 



CSC I 



+ 



l-(^)[-l + 2csc, 



sm I 



2rr 



2 r 



17 15 

— r + ^( 
4 4 



2 1 

cot + 4 CSC + - sin 



2ro 



:ii^ 



For this to be a valid expansion not only must r^/ro << 1 but must also satisfy sin0 >> 
Ts/tq. The tangent to the photon as it travels in the Q = -njl plane of a flat RW spacetime 
is of the form 



k 



con 



urw + cos(0 - Ox - sin(0 - O^w 



(12) 



where urw,Xj 4>RW are respectively unit co- moving time, radial, and azimuthal vectors 
in RW. The significance of the constant angle ^ is that tan^ is the slope of the photon's 
straight line orbit in the co-moving x-y plane, see Fig. 1. When the photon's tangent vector 
is matched across the boundary of the dust hole the following single (exact) constraint results: 



1 + Pb cos (0b - 



n 



sm 



-0 



7o 



(13) 



where r^^cph, and {3^ are evaluated at the photon's entrance/exit point on the boundary of 
the Kottler hole. From Eq. 0131) we obtain the following exact expression for ^ 



sin^ 



where 



\ 



A 
B 



/I 



, (14) 



cos 0b (3b 



To 

rb 



sm 06 Pb — 



sm (pb 
7o 

cos 0b 

7o 



(15) 



The —A choice is made in Eq. f[T^ at the exit point and the Kottler coordinates on the 
boundary are taken as Vb = r2,4>b = 4>2- The +A choice is made at the entrance point and 
the Kottler coordinates are taken as Vb = ri,0b = vr — 0i, (0i is the supplement of the 
entrance azimuthal coordinate, see Fig. 1). 

In what follows we give some of the details necessary to approximately evaluate the 
deflection angle a = ^2 — caused by encountering a condensation in the homogeneous 
dust. The reader not interested in the details can jump to the result called atotai given in 
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Eq. ( 132 p . The calculation is somewhat complicated because the Kottler hole expands as it 
is traversed by the photon. The deflection angle naturally divides into a part that depends 
on the initial size of the hole [astatic given in Eq. fl2^ ] and an additional part caused by the 
extended path required of the photon to exit the expanded hole [agxpand given in Eq. fl?Tl) ]. 
The extended path, described by Ar and A0, is given in Eqs. (128|) and (!30l) . see Fig. 1. We 
will see that the expansion part gives the most significant nonlinear part of the correction 
to the familiar Einstein term AGm/ (?rQ. 



3.1. Approximation Details 

To compute the photon's direction ^ in the dust approximately we assume both ArQ 
and Ts/tq are small (perhaps even of the same order) and expand Eq. (1141) in the two small 
parameters 

and 



3 ro \nj rb 



3 



The result is 



6 

where the coefficients are defined by 



C2(06) = - cot 06 ( i + csc^ 06 ) , (19) 



C3(06) = ^(1-CSC3 06) 



^4(06) = CSC^ 06 



15,„ , . . , , /„ . , 15 . 2 



— (206 - tt) + cot 06 3 - sin 06 - — sin 06 - - sin 06 + - sin 06 + 7 sin 
32 \ ib z o 4 



From Eq. ( IT^ we can conclude that when r^ — )■ 0, sin^ —S = — a/At^Ts exactly with no 
dependence on r6 or 06. This limit is consistent with Eq. (JT8l) . The conclusion is that when 
Tg = there is no A bending. This is an obvious conclusion because the spacetime inside 
and outside of the hole is exactly the same, i.e., no physical difference inside and out exists. 
The only difference is in which coordinates are being used. 
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When < << To we proceed by eliminating ri and r2 using Eq. f[TT|) and then 
expanding 02 about 0i by writing 

02 = 01 + A0. (20) 

This gives us two terms to evaluate 

"static = 6(01 )- ^i(7r - 01 ), (21) 

and 



C'^expand 



+ 



1 



6 



(A0)3 + O (A0)' 



(22) 



The first term can be evaluated immediately using Eqs. (ITSl) and (fT9|) giving the A indepen- 
dent expression 

2 



static 



-2( — )cos<,, 



2 + sin' 



+ 



2rr 



f (201-vr) 



cos 01 ( 4 



15 



sin 01 + 2 sin^ 0i + 7 sin^ 0i + 2 sin^ 0i 



O 



2ro 



(23) 



By overlooking the expansion term one would obviously conclude that there is no A bending. 
To evaluate Oexpand, the expansion's contribution to bending, requires that we compute A0 
(or equivalently Ar) caused by the expansion of the Kottler hole that took place during the 
time it took the photon to transit the hole. In Fig. 2 we indicate how we compute Ar. We 
start by giving the entrance radius ri and look for the common solution to the boundary 
expansion ri,{t) and the photon's radial coordinate rp(t), i.e., we put 



/rr2 



dn 

dt 



dr = c 



dvp 
dt 



-1 



dr. 



(24) 



We rewrite the time it takes the photon to cross the hole as the sum of the time it takes 
to cross from ri on the left to ri on the right plus the extra time it takes to go from ri on 
the right to r2 = ri + Ar. We then move this last time difference to the left hand side and 
obtain the following equation to solve 



ri 



dn 



-1 



/ dTp 
V dt 



-1 



dr = 2c 



ri 



dVp 
dt 



dr = 2c 



7r/2 



dt 



-1 



(25) 



The right hand side RHS is evaluated approximately using Eqs. fIlOp and (fTTj) to obtain 



RHS 



2ro< cot 01 + 



2ro 



cot I 



2 CSC I 



2 log tan ■ 



+ 



Ar| 
18 



cot • 



1 + 2 csc^ 01 



O 



2ro 



+ Ar, 



(26) 
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We call the two terms on the left hand side of Eq. fl25|) LHSb and LHSp and evaluate LHSb 
by expanding in Ar 



LHSh 



1 d f^l 



ri 



2drb \l3bJr 



Id"" 



ri 



6drl V A 



7.^ 



(27) 



where /3b and 7;, as functions of r^, are defined in Eq. ([8]). Equation ffTTj) can be used to 
convert Ar into A0 resulting in 

2" 



Ar 



ro — cos (pi CSC (pi 

1 
2 



2ro 



1+ l_4csc0i 



2rn 



A0 + 



esc' 01 ( 2 - sin" 01 ) + C ( ^ 



(A0)^ + 



1 

+ 6 



cot 01 CSC 0i(l — 6 CSC 0i) + O ( — 

2ro 



(A0)3 + O(A0)n. 



(2^ 



The second term on the left hand side of Eq. (1251) can be evaluated by using 0p(t) from 
Eq. ffTOj) rather than rp{t) (just as was done with RHS) and gives 



LHSp = ro csc^ 0i 



l + O 



+ Arl 



( A0) - cot 01 



(A0)2 + O(A0)H.(29) 



3.2. The Resulting Deflection 



Combining Eqs. (1261) . fl271) and fl29|) in Eq. fl25|) we obtain the change that occurs in the 
exiting value of 0, i.e., A0 = 02 — 0i, caused by the expansion of the hole's boundary as 
the photon traverses 



A0 



— 2/3i sin 01 + 



-/3i I 2 + - sin" (/.I 



3 cos 01 sin 01 
6 sin^ 



2 log \ tan — )■ tan (pi sm (/)i 



+ 



-i/3i Ar2sin0i + O(^ + Ar2 
y \ ro 



(30) 



where f3i is the expansion velocity (f/c) of the dust as seen by observers (who are stationary 
relative to the condensed mass) at the time the photon enters the Kottler hole, see Fig. 
1. Inserting this into Eq. f l22|) we have the additional deflection angle Oexpand caused by the 
extended trajectory of the photon in the Kottler void 



expand 



2rr 



COS 01 



6 sin 01 — 12 cos 0i sin 0i - 



'Ar2 



H sm^ 01 

3 ro 



Ar 



8 20 



■ sm 01 



+ 



+0 



6 cos 01 1^4 sin 0i — sin 0i + 2 sin" 

5/2 



11 sin^ 



12 log { tan y 



Combining Eqs. ( l23l) and (pTl) we obtain the total bending angle atotai caused by a photon 
entering and exiting a Kottler condensation 



«total 



2r, 



cos 01 





2ro, 



2 ~ ~ ~ /Arg 

-4 cos 01 — 12 cos 01 sin 0i ' 



rs ,7 .2/^8 20 2 
— sm'' 01 + Aro - - — sm 



~ / 33 

— (201 - tt) + cos 01 ( 4 H sin 0i - 4 sin^ 0i + 19 sin^ 0i - 64 sin^ 0i 

4 \ 2 



— 12 log <( tan y J> sin^ 0i 



O 



5/2 



The reader should observe that a negative contribution to the bending angle is towards the 
lens and a positive is away. Also recall that these approximate expressions were derived 
assuming sin0i >> r^/ro. Our deflection angle accounts for the finite time (equivalently 
range) that gravity has to act on the passing photon as well as aberration effects caused 
by switching between moving observers. A finite range is equivalent to a shielding of the 
perturbation's mass by the homogeneous distribution of its neighbors, i.e., beyond of 
Eq. ([7]) the effect of the neighbors completely suppress effects of the inhomogeneity. The 
deflection angle Ototai appropriately vanishes in the limit 0i — )■ tt/2, i.e., when the photon 
only grazes the condensation, and for small 0i the lowest order term in the bending angle 
approaches the Einstein value AGm/c^rQ as expected. For an arbitrary impact 0i, however, 
the linear term in atotai is 

ttlincar = "4 ^^^^ 0i, (33) 

and, in some weak lensing circumstances, predicts po tentially detectable differences from the 
Einstein value. In standard lensing calculations [see iBourassa fc Kantowskil (119751 )] there is 
no attempt to make the deflector's gravity fleld part of the cosmology's gravity fleld as GR 
really requires. Deflector masses are simply taken as additions to the cosmology's mean mass 
density and consequently have 'oo' range. Swiss cheese, the only known and relevant GR 
solution, makes the deflector mass a contributor to the cosmology's mean density and as a 
consequence, the gravity fleld of the deflector is limited in range. This limited range is seen 
to be important when the the impact angle is above a tenth of a radian (Figure 4). 



Another somewhat surprising result is that the lowest order correction to the Einstein 
value, other than the flnite time effect represented by the dependence on 0i, is the depen- 



- 11 - 



dence on the expansion rate, i.e., the square root term in Eqs. fl3T]) and fl32l) [see Eq. f|T6l) ]. 
We can interpret the source of this term as the extra time (or equivalently distance) the 
Schwarzschild mass has to act on the passing photon. The Kottler hole expands in size as 
the photon traverses, and since the cosmological constant contributes to the Hubbl e expan- 



sion it contributes to the extra time. Others have also argued that A effects a, e.g.. lSereno 



( I2OO9I ) finds a A contribution to small angle bendin g of or d er (r.s/rnlAr n whic h we do find 
even if of opposite sign and differing amount, llshaki (120081 ): llshak et al.l ( l2008l ) find a term 
of order Aror^ ~ Argcsc^i which we do not. The most important A correction we find, 
i.e., the square root term in Eq. f l32|) . seems to have gone undetected by others because of 
the approximations they used. In the next section we estimate just how important these 
corrections to the Einstein result can be. 



4. Discussion 

In Fig. 3 we have plotted three sets of bending angles for three deflecting masses rang- 
ing from a large galaxy mass to a rich cluster mass, respectively IO^^Mq (lower in red), 
10^'^ M0 (middle in green), and IO^^Mq (top in blue) all at redshift z = 1. Note that red- 
shift z plays a part because redshift influences the entrance/exit size of the Kottler hole, 
see Eqs. and ([7]). The cosmological parameters we used are Qm = 0.3, Q\ = 0.7, and 
Hq = 70kms~^Mpc"^. For each mass we have plotted four bending angles in arcseconds as 
functions of 0i (the supplement of the azimuthal impact angle). The thick lines are |atotai| of 
Eq. ( 132|) . the short dashed lines are | ecstatic | of Eq. (123|) . the dashed lines are ttexpand of Eq. ( 13T]) . 
and the thin solid lines are the Einstein values 2rs/rQ. All deflection angles are negative (at- 
tractive), i.e., towards the deflector, except ttexpand which is away from the deflector. Because 
of the log-log scale it was necessary to plot absolute magnitudes i.e., |atotai| = |«static|— ^expand- 
The reader can easily see (to the accuracy of the plot) that if Oexpand is neglected the de- 
flection angle follows 2rs/ro out to ~ 30°, however, when ttexpand is included the deflection 
angle follows 2rs/rQ only out to ~ 10°. This observation is clearly independent of the masses 
shown and in fact is quite independent of the deflectors' redshifts. The fractional differ- 
ence in atotai and the Einstein value plotted in Fig. 4 is independent of the deflector's mass 
(10^^ — IO^^Mq) and redshift {0 < z < 2) for the range of 0i plotted. Noticeable redshift 
dependent differences would begin to appear for the three masses only below 0i ~ 2°. From 
Fig. 4 we can conclude that for angles above ~ 4° the fractional differences of ttexpand and 
the Einstein values are greater than 1% and above ~ 40° the differences are above 100%. 

In Tabled] we use our corrected bending angle Eq. fl32|) to estimate corrections in bending 
angles for strong and weak lensing by clusters and elliptical galaxies. We look at the following 
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cases: the large image separation cluster lens Abell 1689 at z = 0.18, the high redshift cluster 
RDCS 1252—2927 at 2; = 1.24, and a typical z = 0.5 elliptical galaxy. In A1689, we calculate 
the bending angle corrections for the largest arc separation of 45" for strong lensing and 
weak lensing measurements at 10' away from center (Umetsu & Broadhurst 2008) by using 
the mass profile of recent X-ray measurements (Peng et al. 2009). We also calculate the 
correction in a high redshift cluster RDCS 1252—2927, where the weak lensing signals have 
been detected out to 3' (Lombardi et al. 2005). For lensing by galaxies, we choose a typical 
elliptical galaxy aX z = 0.5 and use the mass profile and weak lensing detections in Gavazzi 
et al. (2007). In general, we find the corrections in the bending angles for strong lensing are 
quite small, e.g., the largest correction (1 — cos^^i) is just 0.07% for the largest separated 
arcs in A1689. However, for weak lensing, the correction can reach 2% for the weak lensing 
signals detected in the outermost regions of the cluster in A1689, and the correction can 
reach 5% for the z = 1.24 cluster RDCS 1252—2927. For the weak lensing signals detected 
using an ensemble of elliptical galaxies (Gavazzi et al. 2007), the correction is 4% for the 
outermost bin. A correction of this amount will present an additional challenge for using 
weak lensing as a tool for precision cosmology. For the corrections involving the A term, 
the largest is 0.02% for weak lensing in high redshift clusters which is not detectable in 
current observations. We expect our model to be relevant for weak lensing induced by the 
large scale structure including weak lensing of the cosmic microwave background, where even 
larger volumes are involved. We expect a large correction due to the 1 — cos^ term and a 
presumably detectable correction involving A. 

The corrections we give for strong lensing are negligible because (pi is small and only 
a small fraction of the inhomogeneous mass appears inside the Einstein ring. For these 
cases, our corrections may not be accurate because the effective lensing mass is not spher- 
ically distributed as it is in our model. Non-linear corrections are conceivably sensitive to 
the difference in cylindrical and spherical symmetry. More realistic models are needed to 
fully constrain corrections for strong lensing. In general, the applicability of the corrected 
deflection angle atotai in Eq. f l32|) is limited to spherical inhomogeneities, the majority of 
whose mass is within the minimum impact of the light ray. This is because we used a fully 
condensed Swiss cheese model, i.e., the homogeneity is represented by a Schwarzschild mass. 
Because we are calculating non-linear corrections one cannot expect Eq. fl5^ to give an ac- 
curate answer by simply including that fraction of the mass within the impact cylinder as 
is normally done in lensing. Consequently more accurate mass profiles in the Swiss cheese 
would be appropriate for the strong lensing examples in Table 1. 

Work on this paper was initiated to correctly quantify the cosmological constant's ef- 
fect on small angle deflections of photons caused by mass inhomogeneities in an otherwise 
homogeneous cosmology. By using an exact solution to GR we established that A's effect is 



- 13 - 



non-linear thus requiring use of a gravity theory beyond Newton's. The model we used, a 
flat Swiss cheese cosmology, also predicts a significant decrease in the deflection angle caused 
by the shielding of an inhomogeneity by its homogeneously distributed neighbors. Shielding 
occurs because the deflectors are contributors to the cosmology's mean density. Standard 
lensing calculations completely overlook shielding because deflectors are treated as additions 
to the mean. 

Perturbations to Ototai would obviously exist if the neighbors generated a shear at the 
site of the deflector. The accuracy of Swiss cheese predictions depends on the scale at which 
inhomogeneous matter follows the background Hubble flow, i.e., on what scale the cosmo- 
logical principle is satisfied. The simple Swiss cheese model used here doesn't allow for 
peculiar motions but does account for the scale of the cosmology effected by an inhomo- 
geneity, i.e., beyond the perturbed spacetime returns to the mean cosmic flow. In the 
neighborhood of the Local Group, where good observational data is available, most galaxies 
follow the Hubble flow with only small deviations (e.g., Karachentsev et al. 2009). A hierar- 
chical Swiss cheese condensation could be used to include shear and peculiar motion but it 
would not only complicate this calculation by introducing several additional parameters, it 
would most certainly obscure the source of the A term in the results. To keep the result as 
simple as possible we did not attempt to estimate the size of these additional perturbations. 

Our results e.g.,Eq. f p2|) are stated in terms of the parameters cpi and tq described in 
Section 3 and Figure 1 and are not necessarily the most convenient ones to use in lensing 
applications, however, they were convenient for the above derivations. To have the incoming 
photon travel parallel to the x-axis one only has to rotate the coordinates clockwise an 
amount given in Eq. (fT^ . 

This work was supported in part by NSF grant AST-0707704 and US DOE Grant 
DE-FG02-07ER41517. B. Chen also thanks the University of Oklahoma Foundation for a 
fellowship. 
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Table 1. Examples of Gravitational Lensing Corrections in ACDM Cosmology 



Name 


Lensing 


redshift 


Mass 




Impact Angle 


01 


1 — cos3 <f)i 


ratiol^ 


ratio2'^ 










(Mpc) 


(arcsec) 


(degrees) 








A1689 


strong 


0.18 


8 X 10^3 


6.6 


45 


1.2 


0.00065 


2.2 X 10-^ 


1.4 


A1689 


weak 


0.18 


10i5 


15.3 


600 


6.8 


0.021 


0.00017 


1.4 


RDCS1252-2927 


weak 


1.24 


1015 


8.0 


180 


11 


0.052 


0.00040 


0.20 


Elliptical Galaxy 


strong 


0.5 


3 X 10" 


0.8 


2 


0.87 


0.00035 


1.7 X 10^'^ 


0.69 


Elliptical Galaxy 


weak 


0.5 


1013 


2.6 


70 


9.6 


0.041 


6.6 X 10-5 


0.69 



^ratiol = 4 tan 01 y + ^ sin^ (f>i [the ratio of the next order correction to the lowest order term, see Eq. I|32|l ] 
^ratio2 = ^= ^ „ - (measures the relative importance of A in the square root term). 
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Fig. 1.- 

A photon travels left to right entering a Kottler hole at r = ri, = tt — 0i and returns to 
the FLRW dust at r = r2, = 02- The photon's orbit has been chosen symmetric in Kottler 
about the point of closest approach r = ro,0 = 7r/2. Due to the cosmological expansion 
Ar = r2 — ri > 0. The slope of the photon's co-moving trajectory in the x-y plane is 
when incoming and ^2 after exiting. The resulting deflection angle as seen by a co-moving 
observers in the FLRW background is ct = ^2 — Ci- 



y 




Figure 1 
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Fig. 2.- 

While a photon travels through a Kottler hole entering at r = ri, = tt — 0i and exiting at 
r = r2, = 02 its radial coordinate varies with time as rp{t) and the boundary of Kottler 
hole continually expands according to rb{t). The exit coordinates differ from the entrance 

values by A0 — (f)2 — (pi and Ar — r2 — ri. 

r 





1 


nit) 


Ar 


i 




rM) ^ 








ri 


^0 


At 





Figure 2 
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Fig. 3.- 

Three sets of the deflection angles a corresponding to three deflector masses lO^^M© (lower 
in red), IO^^Mq (middle in green), and IO^^Mq (top in blue) at redshift z = 1 are shown as 
functions of the azimuthal angle 0i. The thick lines are |atotai| of Eq. (132|) . the short dashed 
lines are | astatic! of Eq. (1231) . the dashed lines are ctexpand of Eq. (I3T]) . and the thin solid lines 
are the Einstein values 2rs/rQ. All deflection angles are towards the deflector except aexpand 
which is away from the deflector, i.e., in this plot jatotail = |«static| — ^expand- 




Figure 3 
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Fig. 4.— 

The fractional difference of the Einstein deflection angle 2rs/rQ and atotai given in Eq. (15^ . 
i.e., {2rs/rQ — | atotai I ) / 1 ^totai I as a function of (pi in radians. For the domain of (pi plotted 
and to the accuracy shown, the fractional error is remarkably independent of the mass of 
the deflector (lO^^ - lO^^M©) and its redshift (0 < ^ < 2). 




02 (radians) 



Figure 4 



